Theory of the Quantum Critical Fluctuations in Cuprates 
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FIG. 1: Current pattern(a) and its abstraction (b) 

in fig. (lb). A classical model with such a pair of Ising 
degrees of freedom is the Ashkin- Teller 8] model: 



H = — ) Ji{(Ji(Jj 



JiiOiTiGjT] , 
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The statistical mechanics of the time-reversal and inversion symmetry breaking order parameter, 
possibly observed in the pseudogap region of the phase diagram of the Cuprates, can be represented 
by the Ashkin- Teller model. We add kinetic energy and dissipation to the model for a quantum 
generalization and show that the correlations are determined by two sets of charges, one interacting 
locally in time and logarithmically in space and the other locally in space and logarithmically in 
time. The quantum critical fluctuations are derived and shown to be of the form postulated in 1989 
to give the marginal fermi-liquid properties. The model solved and the methods devised are likely 
to be of interest also to other quantum phase transitions. 



A mean-field solution of a microscopic model predicts 
that the pseudogapped region [l| of the phase diagram 
of the cuprates (Region II in fig. (1) of Ref. (2)) has 
a spontaneous current pattern depicted in fig. (la) Q. 
Such a state breaks time-reversal, inversion and three of 
the four reflection symmetries of the square lattice while 
preserving translational symmetry. The observed mag- 
netic diffraction in neutron scattering 3] as well as the 
dichroism observed in ARPES Q with circularly polar- 
ized photons is consistent with such a state. 

The transition temperature of this state varies contin- 
uously with hole density x and tends to at x = x c , the 
quantum critical point. Anomalous normal state proper- 
ties are observed in the funnel shaped region (Region I in 
fig.(l) of Ref. (2)) emanating from the QCP. They were 
shown to follow from a phenomenological spectrum ||: 

Imx(q,w,T) ex u/T, for u « T, (1) 
oc constant, for T«w<C u> c , 

where oj c is a cut-off, determinable from experi- 
ments. The corresponding real part of the fluctuations, 
i2ex(q, co,T) oc hx[u) c /ui) for oj/T 3> 1 and oc hx{uj c /T) 
for uj/T <C 1. Since this spectrum has a singularity in the 
limit (T, uS) — ► 0, it specifies the fluctuations near a quan- 
tum critical point. For small q, this spectrum is directly 
observed in Raman scattering. Jjj. The two unusual 
properties of EqQ] lo/T scaling pfl and no dependence 
(or a smooth dependence) on q are sufficient to gener- 
ate a marginal fermi-liquid ||. The observed anomalous 
properties are well explained by the marginal fermi-liquid 
and its predictions for the single-particle spectra have 
been verifed in experiments. 

The purpose of this work is to show that the quantum 
critical spectrum of Eq.Q is the spectrum of fluctua- 
tions of the order parameter which condense to give the 
observed order depicted in fig. (la). 

To do so we consider a quantum generalization of the 
classical model whose solution gives such an order pa- 
rameter. The four possible configurations in each unit- 
cell of fig. (la) can be represented by four vectors with 
arrows representing time-reversal and orientation repre- 
senting the only plane of reflection symmetry, as shown 



(2) 



o~i and T^are Ising spins. Such a Hamiltonian can also be 
derived |2j from the microscopic model which in mean- 
field theory gives the observed order parameter provided 
one neglects amplitude fluctuations. 

The Ashkin- Teller model has a variety of phases de- 
pending on J 4 /J 2 and T/J 2 HHE3- Especially inter- 
esting to us is the region —1 < J4/J2 < in which the 
ordered phase has ferromagnetic order in < o~i > , < t% > 
as well as < OiTi > just as in the order observed. Two 
more properties of this region are of special interest: on 
the critical line, the model is Gaussian (with stiffness de- 
pending on J4/J2) and the specific heat at the transition 
has no divergence. 

With the transformation 

cos(0i) = fa + Ti)/2; sin(^) = fa - n)/2, (3) 
the Ashkin- Teller model becomes 



H 



i,jnn 



2.J2 cos( 



J4C0s2(#i 



) + hi cos A6i. 

(4) 

The restriction that 6^ = 0, 7r/2, 7t, 37r/2 has been written 
in terms of a four-fold anisotropy term oc ft, 4 . Monte- 
Carlo calculations show [ll| that the classical phase di- 
agram of model J5J and of are the same. This is 
consistent with the fact that h\ is irrelevant above the 
critical point in the xy model; hi however is relevant 
below the transition enforcing an Ising state |l2j . 
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For a quantum theory for the model, two additional 
features are included, the kinetic energy and the dissipa- 
tion. For the latter, we consider physics of the Caldeira- 
Leggett type 0], which can be derived by integrating 
over the fermions in the microscopic model which cou- 
ples the orbital-currents of (8i — 8j) to the currents of 
the fermions ^4| • For simplicity of presentation, we start 
with J4 = 0. We will show that no essential modification 
arises at least for | J4I/J2 < 1/2, which is in the interest- 
ing regime. We will also take /14 = and show that in 
the fluctuation regime, which is our primary interest, a 
finite I14 does not change the results. With these simpli- 
fications, 



S = S Q + J dr2J 2 cos (fii (r) - dj (r)) , (5) 

7*?T(*)' + 

J dT j dT , 



5', 



(r-rO 



where % (r) = 0i (r) — 8j (r). 

We use the Villain representation of the periodic func- 
tions and introduce integer link variables = rrii X x + 
rriiyy and discretize time in steps At to get 



(ij) ^ 



— Bi V ) + S, 



diss 



(6) 



After a Fourier transform and integration over the 9i Tl 
variables, the action for the vector field m is obtained: 

s = 7^E( 27r ) 2G ( k <' w »)k a4 i k < xm <«i 2 



J1 



a k l \ „2 / a _s2 , , _ ,2 



Aire \ujr, 



a (At) \uj n mi n \ 



(7) 



Here k; are discretized momenta, k; = (2,Trl x /a)x 
+ (2irl y /a)y, uj n are the Matsubara frequencies, oj n = 
27TO//3, and J = 2J 2 Ar/a , J t = 1/EaAr. The prop- 
agator G (k;, uj n ) is given by 



G Qci,w n ) = ( Jckf + Jttol/c + a \u n \ kf/An)' 



(8) 



Quantum dynamics introduces sources and sinks in the 
vector field. Therefore this action includes besides the 
vortices, (curl m), the time-derivatives of m. Sources 
and sinks due to quantum dissipation suggest that the 
time-derivative must include a field with a divergence. 
We define two orthogonal charges, Pt,(r, t) and po( r ;t) 
(in the continuum limit) through 



p v (r,t)z = V x m(r,t), 



(9) 



Po (r,t) = (l/c)dm(r,t)/dt-T 



(10) 



where the velocity c = a/ At. The time derivative of m 
appearing in the action then is 



2 1 \2 



w n 1 |2 . 2 I |2 



(11) 



In the theory of the xy model, the charges p v are sources 
of "vortices" where the velocity field is azimuthal with a 
strength falling off as 1/r. po are sources of the orthog- 
onal radial field ^(|. From IjlOfl. it is seen that po's are 
events in time where the divergence of the field m(r, t) 
changes. The importance of po lie in their logarithmic 
interaction in time. 

In terms of p v ,Po, the action, in the continuum limit, 
neatly splits into three parts: S = S v + Sq + S' 



So 

D 



J 

\p v {kiu n )Y 



|p (ki^«)| 2 



— y 
—y 

x \p (kiuj n )\ 2 



(12) 



Aire 



a 2 kf 
16tt 2 



The interesting part about this decomposition [15| is 
that S v is the Coulomb gas representation of the xy 
model with charges interacting logarithmically in space 
but locally in time and Sq is the Coulomb gas repre- 
sentation of the Kondo problem with charges interacting 
logarithmically in time but locally in space. S has non- 
singular interactions and is unimportant compared to S v 
and So H3- Being orthogonal the charges p v and p 
are uncoupled; the action is a product of the action over 
configurations of p v and of p . Any physical correlations 
are determined by correlations of both charges. Both the 
(PvPv) (k, u>) and the (poPo) (k, oj) correlations are well 
understood. For any finite J, the field p v is confined 
in the limit T — > 0; no free vortices exist. The phase 
transition can come about only due to free po's, due to 
a tuning of dissipation parameter a. We therefore first 
remind ourselves of the correlation function of the po's. 

Let us introduce a core-energy A for the po's just as 
is done to control the fugacity of p v 's, the vortices. Next 
consider how the renormalization of a and A proceeds. 
Including the core-energy, the action Sq is 



So = y 



^E 



47T \uj n \ 



\P0i (Wn)T 



^tA \p 0l (t) 



(13) 

This is identical to the Coulomb gas representation 
[lij l of the Kondo problem and a quantum dissipation 
problem 0|. The resulting RG equations are 



dA 



(l-a)A 



(14) 
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da ~n 



where A = Ar c and r c is a short time cutoff. The critical 
point of interest is at a c = 1 |2l|, where A scales to 0; 
for a < 1 the charges freely proliferate as "screening" 
due to A becomes effective, a > 1 represents the or- 
dered or confined region in which the anisotropy field 
is strongly relevant. We are interested here only in the 
region a < 1. Well in the quantum critical region, the 
(singular part of the) propagator for po is 



(Po (w n ) Pa (-^n)) = (l/4?r \u n \ 



(15) 



The crossover to the quantum-disordered or screened 
state is given when w x is of the order of the inverse of the 
characteristic screening time, which may be estimated 
similarly to Kosterlitz's estimate of the screening 
length in the xy-problcm. Thus 



r c 1 cxp (—b/y/l — a) , 



(16) 



where b is a numerical constant of 0(1). At finite tem- 
peratures and low frequencies, the crossover temperature 
may be taken as given by Eq. Ijltil ) with u x replaced by 
T 

Finally, we come to the correlation function of interest, 
that of the order parameter: 



(7 ■ — le lBi »p- 



i,) = 111 



fd[0], 



—s 



(17) 



To compute this we will employ a procedure similar to 
the one used by Jose et al. |l3| for the 2d xy model. 
Consider a path from ip to jv, and split it into two: 
ifi — > jp and jp — > jv. All paths should give the same 
answer so we have chosen the one most convenient. A 
vector field lf ifl is defined which lives on the sites of the 
lattice and whose components are 1 if the path crosses 
site ip and zero otherwise. To capture the second part of 
the path, we define a scalar field rjf^, which lives on the 
links in time, and is nonzero only for paths on the link. 
Including the fields described above and integrating over 
the 9's, we get 

S = S + -^]T(27r) JG(k l ,uj n )[(uj n 7 1 */2)(tk 1 - min ) 



In 



+ (k/ • lf*/2) ■ m in) + c.c. 
+ L„77 /2 + k ; -lfV2| 2 ]. 



(18) 



The last term in the sum gives the spin wave contribu- 
tion to the correlation function. Clearly the spin wave 
contribution is finite and this is indeed the standard re- 
sult that spin wave fluctuations do not disorder the state 
in three dimensions. The second term in the sum can be 
written in terms of the linear coupling of (k x if) to p v 's, 



the vortices. The vortices are confined for T — > for the 
case of a finite J assumed by us and are therefore not 
interesting. The first term can be written in terms of a 
linear coupling between 77° 's and po' s - This is the inter- 
esting term because the divergence in frequency induced 
by the dissipation as a is tuned leads to a proliferation 
of po's and thus to disorder. The correlation function 
near the quantum critical point are thus determined en- 
tirely by the dynamics of po's. The correlation function 
is proportional to, 



C(y-y',t-t') = exp(-F), (19) 
F = -^TV f dk^ (p a K) po (-«„)) x 

G 2 (k, w„) (1 - cos (k • (r - r') + u n (r - r'))) . 

Deep in the quantum critical regime the correlations 
function (po {ui n ) po (— w„)) is given by Ea. H15|) . The 
summand over n in F has a leading l/|w| part. It 
is easy to see that for any finite |r — r'| the sum over 
71 is divergent. Thus, for any finite spatial separation, 
the correlation function C(y — y',t) is identically zero 
for r larger than the crossover scale r c . On the other 
hand, for |r — r'| = 0, the correlation function is given 
by C(0, t-t')= exp(-F(0, t - r')), where 

win >\ 9 1 - cos (tJn (t ~ t')) 

F(0,t-t ) = -2ttT > 1 1 log(|u; n |T c ). 

(20) 

Such correlation functions have been calculated in 
other contexts. In particular Ghaemi et al. 22], pro- 
vide the spectral function of the correlation to be 



Jmx(q, ui) = t c tanh(w/2T). 



(21) 



where Tc = cor c , cq is a constant of O(l). The phe- 
nomenological spectrum of Eq.QJ is thus derived. 

We now summarize the calculations for the effect of 
the anisotropy field hi. As in Ref . ( fl^| ) . we supplement 
the Action of Eq. 112|l by introducing a new field pi 



hi cos(4fli) 



(22) 



where j/4 = /14/2. Introducing the m fields as before 
and performing the 6 integral leads to additional terms 
in the action which couple pi linearly to po.i- This lin- 
ear coupling can be eliminated with a renormalization 
of the coupling constant a to 2a. This does not change 
the critical behavior in the quantum critical regime; the 
correlation remains of the form (|21|l . For a > a c , hi is 
relevant just as in the classical problem and enforces an 
Ising state. 
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The coupling J4 is easier to treat. For —0.5 < J4/J2 < 
0, the bare potential continues to have the absolute min- 
ima at (9i — 6a) = 0. Therefore the Villain representa- 
tion of periodic functions can be made as above with an 
altered coefficient. This affects the cut-offs in the solu- 
tion but not the form of the correlation function below 
the cut-offs. The theory is valid only in this range; for 
smaller J4, an Ising transition to a different symmetry is 
expected as in the classical model. 

We summarize the principal results: We have consid- 
ered a statistical mechanical model which has the sym- 
metries and the degeneracies of the observed phase in 
underdoped cuprates. The specific heat at the classical 
transition in this model is a smooth function of temper- 
ature as in experiments. We have generalized the model 
by including inertial as well as dissipative dynamics. On 
tuning dissipation, the model has a phase transition at 
T = 0. The critical fluctuations of the model are de- 
termined by the dynamics of charges with interactions 
which are spatially local but logarithmic in time. The 
correlation function in the quantum critical regime calcu- 
lated for the model has the form postulated phenomeno- 
logically to understand the properties of cuprates in the 
quantum critical regime, i.e they display uj/T scaling and 
spatial locality. 

More generally, our results are applicable also to the 
dissipative quantum xy model and hence relevant to criti- 
cal properties of Josephson Junction arrays in two dimen- 
sions. The Ashkin- Teller model is a staggered 8-vertex 
model. We expect that our method has application to 
quantum versions of 6 and 8-vertex models generally to 
which many physical problems of interest correspond. 

We wish to thank T. Giamarchi, B.I. Halperin, A. She- 
hter, A. Sudbo, and A. Vishwanath for useful discussions. 
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